Spin-orbit Coupled Bose-Einstein Condensate under Rotation 
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We examine the combined effects of Rashba spin-orbit (SO) coupling and rotation on trapped 
spinor Bose-Einstein condensates (BECs). Nature of single particle states is thoroughly examined 
in the Landau level basis and is shown to support the formation of half-quantum vortex. In the 
presence of weak s-wave interactions, the ground state at strong SO coupling develops ring-like 
structures with domains whose number shows step behavior with increasing rotation. For fast 
rotation case, the vortex pattern favors triangular lattice, accompanied by the density depletion in 
the central region and weakened Skyrmionic character as the SO coupling is enhanced. Giant vortex 
formation is facilitated when SO coupling and rotation are both strong. 
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Introduction.- Bose-Einstein condensates (BECs) are 
a fascinating testing ground for many interesting con- 
densed matter physics phenomena such as the formation 
of vortices by rotation. For a single component conden- 
sate, if the rotation is fast enough, the vortices form a tri- 
angular (Abrikosov) lattice, indicating the entry of BECs 
into a quantum Hall regime where the lowest Landau 
level (LLL) physics dominates [1]. In the case of spinor 
condensates, the competition between inter- and intra- 
component interactions was found to yield even richer 
vortex patterns such as square lattice, double-core lat- 
tice, and interwoven "serpentine" vortex sheet [2], or 
even giant Skyrmions [3] if the symmetry between the 
two components is broken. 

A recent focus in cold atom physics is on the influence 
of "engineered" Rashba-type or Dresselhaus-type spin- 
orbit (SO) coupling among spinor BECs induced by the 
so-called "synthetic non-Abelian gauge fields" [4]. The 
SO effects in cold atoms have received increasing atten- 
tion theoretically [5-7] and its first experimental real- 
ization was recently reported [8]. Non-trivial new struc- 
tures, such as stripe phase and half-quantum vortex, have 
been found or predicted, enriching the phase diagram of 
spinor BEC system [6, 9]. 

All existing studies of SO effect in spinor condensates 
refer to non-rotating case, and the aim of this paper is 
to study the combined effects of Rashba SO coupling 
and rotation on spinor BECs, focusing on spin- 1/2 case 
as the prototype. Even the single particle aspect of the 
Rashba BEC Hamiltonian in the presence of rotation is 
not well understood theoretically, and we address these 
issues first. We use the Landau level (LL) expansion to 
examine the ground state of the single particle system 
by exact diagonalization, then move on to discuss the 
influences of weak 5- wave collisional interactions. New 
features which are absent from pure Rashba, or pure 
rotational BECs, are observed. 



tide Hamiltonian describing the trapped Rashba SO 
coupled system under rotation Hq = /^'''(p^ + 2/^p • 
o-)^/(2m) - QzLz + /^V(r)^, where p = -ihV 
with m being the particle mass, a are the 2x2 Pauli 
matrices, and hz denotes the strength of SO coupling. 
Note that we make implicit the 2x2 unit matrix for 
p. ^ = (?/;^,?/;^)^ describes the two-component wave 
function. Qz stands for the rotational frequency and 
Lz = /^^(r X p)^^ is the canonical angular momen- 
tum, while the non-canonical part oc kt x cr/m can be 
compensated for an external spatially dependent Rabi 
coupling between the two components [10]. The trap- 
ping potential V^(r) = muj'^{x^ + y^)/2 is introduced to 
keep the boundary condition realistic. We assume that 
the trapping frequency in the z direction is much higher 
than those of transverse directions, so we may view our 
condensate as an effective two-dimensional system with 
the degree of freedom in the z direction frozen out. With 
the proliferation of parameters in the problem, it is con- 
venient to turn all quantities dimensionless by introduc- 
ing a length scale ^Jll|muJ and an energy scale huj. The 
dimensionless Hamiltonian thus becomes 
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zV + 70- - 6i X r)^| 



+76 l^^^z X r • o-)^ + -(1 - \?) Ir^"^^^. (1) 

Two dimensionless parameters 7 = n/yh/muo and h = 
Vtz/^ are introduced. Note that the second term of 
Eq. (1) is induced by the interplay between rotation and 
SO effects. We require 6 < 1 in order to prevent the 
complete cancelation of the trapping potential by the 
centrifugal force. 

The solution of the single particle Hamiltonian (1) 
is conveniently obtained from its second-quantized form 
which reads 

7t- 



HQ = {l^h)d:a^{l-h)ffc^i^ 
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Single Particle States.- We start with the single par- where operators are introduced as a 



dz. o) 




FIG. 1. (Color online) Single particle ground state of rotating 
Rashba-BEC Hamiltonian. (a) Average angular momentum 
{c^ c) and (b) average LL index {o/ a) for -0^ in (6, 7) parameter 
space. Cubic roots ((c^c)"^/^ and {d)aY^^) are plotted for 
improved view, (c) to (e) are radial density distributions of 
V^^ (black dotted) and -0^ (red solid) for the same large 7 and 
increasing h. The numbers with arrows denote the winding 
numbers of wave functions, (f) Densities of both components 
for the same b and increasing 7. The digits in the legend label 
the winding numbers. 



^z — dz^ c = ^z -\- dz, and c^ = ^z — dz in complex 
coordinates z = x -\- iy^ z = x — iy^ dz = {dx — idy)/2^ 
and dz = {dx -^idy)/2. The two pairs of operators satisfy 
the bosonic commutation relations, [a^a)] = [c^c^] = 1. 
We focus on positive b and 7 without losing generality. 
Without SO coupling, each Landau level state would 



be an eigenstate satisfying d^d(j)n,' 



and 



c^c(t)n,m = '^<j>n,m with n and m being the Landau level 
index and angular momentum number. The eigenenergy 
is £0 = (1 + ^)^ + (1 — b)m^ i.e., (1 + b) and (1 — b) 
are the energy spacing between LLs and angular mo- 
mentum states in each LL, respectively. The ground 
state is achieved for n = m = 0. For nonzero 7, on 
the other hand, due to the specific form of Rashba term 
which couples neighboring LLs, one can use the ansatz 
that the eigenstate be formed as the linear combination 

^ = ^m>ui^rn(l)m-u+l,m, yni(t>m-v,niY - By Substitut- 
ing ^ into Hq^! = e^f with e being the eigenenergy, we 
find that the coefficients Xm^ym obey 

[(l + 6)(m-i/+l) + (l-b)m]xm 

+ i-f{^/m-v^lym - VmTT^m+i) = exm, 
[{l^b){m-v)^{l-b)m]ym 



The one-dimensional coupled linear equation set can be 
solved by numerical exact diagonalization with a reason- 
ably large cutoff in m for each v. Note that v is an inte- 
ger number defining the phase windings of '^^ and ijj^ as 

V — 1 and V, respectively, as shown in Fig. 1(c)- 1(f). The 
difference by 1 indicates the formation of half-quantum 
vortex [11]. When either 6 or 7 increases, the value of 

V in the ground state increases correspondingly. In the 
limit 6 = 1, each LL becomes infinitely degenerate, and 
we find by numerical approach the divergence of v due 
to the disappearance of trapping potential. 

In Fig. l(a)-l(b) we show the average angular momen- 
tum {c^c) and average LL index {d)d) for ij)^ (ij)^ shows 
similar behavior). The visible steps in figures are strongly 
related to the discreteness of LLs. For fixed 7, increas- 
ing b results in higher angular momentum and lower LL 
index due to the decrease of energy spacing between an- 
gular momentum states within each LL (1 — b) and the 
increase of energy spacing between LLs (1 + 6). It also 
results in the more ring-like structure of the ground state 
density with increasing radius at larger b. By calculat- 
ing the local spin density vector n = ^"^ cr^ , i.e., rix = 

we also find that the spin texture at 6 = corresponds to 
the Skyrmion configuration which can be characterized 
by the Skyrmion density 

1 n • {dxn X dyii) 



Ps{^) 
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(4) 



As b increases, we observe the growing overlap between 
densities of two components, resulting in the decrease 
of ps due to n^ ^ 0. Note that the Skyrmion density 
is strictly zero for any planar spin texture. Along with 
the density depletion in the trap center, the ground state 
exhibits spin spirals in the annular region. 

When SO coupling is absent, b ^ 1~ limit is associ- 
ated with the quantum Hall regime. For small 7, we may 
view the SO coupling as a perturbation and the LLL 
picture still applies. However, when 7 increases to be 
comparable to (1 — 6), higher angular momentum states 
contribute more and more to the ground state and the 
ring-like nature as well as a large winding number u be- 
comes dominant as seen in Fig. 1(f). For the value of b we 
choose, when 7 is comparable to the LL energy spacing 
(1 + 6), the radius becomes enormously large. 

Interaction Effects.- We may now include interactions 
in the manner suitable for spin- 1/2 condensate, 

Hi = ] [[gdPt + Pi? - 9s{Pt - Pi)% (5) 



• iJ{^/rnXm-l - Vm-u^lXm) = eVn 



(3) 



where gc = {g^ + ^n)^ ^^^ 9s = {gn - gn)^ with 
g^^ and g^i being the intra- and inter-component inter- 
actions, respectively. It is assumed that g^^ = g^ and 
gtl — 9 it' ^ i^ ^^^ particle number, thus pi = l^/^ip 
{i =t, D satisfy J ^^ Pi = 1. The full Hamiltonian be- 
comes H = Ho-\-Hi. The ground state is numerically ob- 
tained using the imaginary time evolution method. Now 



we also consider two specific cases as in the single particle 
analysis above. 
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FIG. 2. (Color online) Ground state evolution with b when 
7 = 10 and {gc, Qs) = (9,-1). (a) and (b) are plots of the 
density iV^tT ^^^ phase arg('0^) for b — 0.0, while (c) and (d) 
are for b = 0.02. The long arrows in (b) and (d) represent 
the directions of k vectors, (e)-(i): Up-component densities 
|^^|2 for b = 0.08, 0.12, 0.16, 0.20, and 0.50, respectively, (j) 
Dependence of domain number, total density difference J \nz\ 
(xlO), and angular momentum of ijj^ on b. ipi shares similar 
evolution behavior. 

We start with strong SO coupling case. When there is 
no rotation, the ground state is a stripe (ST) {gs > 0) or 
a plane wave (PW) {gs < 0) as was found recently [6]. In 
our trapped system, we choose weak interactions so as to 
recover the same stripe and PW phases. We first consider 
^s < case and discuss ^5 > later. The PW phase is 
composed of a single k state as shown in Fig. 2(a)-2(b). 
As b increases gradually above zero, the — k state also 
begins to contribute, and each component of the conden- 
sate is divided into two parts with opposite momenta, 
indicated by long arrows in Fig. 2(d). Vortices emerge 
due to the counter-propagation of plane waves. Loca- 
tions of these vortices also coincide with the low den- 
sity regions shown in Fig. 2(c). We notice the similarity 
of this phase to the metastable state found in Ref. [6]. 
As we keep on increasing 6, more k states contribute, 
and a single domain of vortex line (Fig. 2(c)) turns into 
several symmetrically placed domains with different k's 
(Fig. 2(e)-2(h)), while vortices nucleate at the center of 
the trap. The number of different k domains increases 
as3^4^5^6^--- until each component of the 
condensate becomes indistinguishable from a concentric 
ring (giant vortex) at strong enough rotation as shown 
in Fig. 2(i). The ring then expands in radius as b further 



increases. The integer number of domains shows step 
behavior with rotation as shown in Fig. 2(j). Addition- 
ally, the angular momentum of the condensate increases 
monotonously with 6, showing slight jumps (indicated 
by short arrows in Fig. 2(j)) where the domain number 
changes. We also observe the enhanced overlap between 
densities of two components as b increases. We define 
J \nz\ as the total density difference to quantify the over- 
lap, where the absolute value is chosen to avoid cance- 
lation among opposite signs of n^ upon the integration. 
From Fig. 2(j) we observe the obvious overall trend of 
decrease, where the minor turns also coincide with the 
transition points of the domain number. 




FIG. 3. (Color online) Ground state evolution with 7 for 
large rotation b = 0.998 and (^c, gs) = (9,-1). From left 
to right: density of up-component I'^tT? magnitude of its 
Fourier transformation | Jexp[— ik • rJlV^^I^I, total density of 
two components (iV^tl^ + IV^-^I^)? ^^^ ^^e Skyrmion density 
Ps. (a)-(e): 7=0.0, 0.002, 0.1, 0.2, and 0.3, while the total 
Skyrmion number, defined as the spatial integral of Eq. (4), 
is 40.9143, 19.2982, 12.8454, 8.1136, and 3.7456, respectively 

So far, we have discussed ^^ < case. For ^^ > 0, 
the ground state is a stripe phase corresponding to the 
standing wave formed by two counter-propagating plane 
waves when b = 0. As we increase rotation, we also 
observe similar evolution behavior as already shown in 
Fig. 2(c)-2(i). This implies that for weak interactions, 
as b increases, the difference caused by the sign of gs is 
smeared out, giving way to the effects of large 7 and b. 

When there is no SO coupling, the ground state in the 
so-called quantum Hah regime {b ^ 1~) is largely com- 
posed of the LLLs. Previous theories for the symmetrical 



case {g^^ = g^^i) found an interlaced vortex array for each 
component, while the detailed lattice structures depend 
on the value of gs [2]. Starting from the known limit, 
7 = 0, we hereby examine the influence of finite 7. For 
consistency, we choose the same values of interactions as 
in the strong SO coupling case and discuss the infiuence 
of gs later. Square vortex lattice {gs ^ 0~) is shown in 
Fig. 3(a), along with the corresponding Skyrmion den- 
sity Ps(r). Interestingly, 7 = lattice of displaced vortex 
cores of i/j^ and i/j^ already exhibits a regular pattern of 
Skyrmion density, qualifying it as a Skyrmion lattice [12]. 
We emphasize that the same Skyrmion charge is obtained 
for both spin- up and spin-down vortex cores. 

From top to bottom of Fig. 3 we can find that, with 
increasing 7, the vortex lattice would be distorted first, 
and then densities deplete away from the trap center, 
forming the ring-like structure as in the single particle 
case. In the annular region of Fig. 3(d)-3(e) the vortices 
prefer to form triangular lattice pattern, even though we 
started with the square lattice for 7 = 0. Inside the 
trap the giant vortex is formed. During the transition 
from square to triangular lattice, we observe that the 
Skyrmion distribution shrinks first, and then expands 
as b increases as shown in Fig. 3. To characterize the 
trend of the Skymrion density becoming more squeezed, 
we define the total Skyrmion number as J ps which de- 
creases correspondingly (see the caption of Fig. 3). This 
fact results from the increased overlap between up- and 
down-component densities. When the vortices overlap 
completely, Skyrmions are also destroyed. 

Ground state of pure rotating spin- 1/2 BECs may also 
be interlaced triangular vortex lattice when gs is close 
to —gc {g^i is small compared with ^^^), or double-core 
lattice when gs = 0, or even vortex sheet when gs > [2]. 
For weak interactions, as we increase SO coupling, we 
observe similar evolution pattern as shown in Fig. 3 
for all these vortex structures. Ring structure with a 
giant vortex about the trap center would appear, while 
triangular vortex lattices is favored in the annulus, 
and the overlap between densities would be enhanced. 
As in the case of strong SO coupling discussed above, 
the details of gs play no qualitative role for large 7 and b. 

Summary.- We examined several prominent features of 
the rotating trapped Rashba BECs by analytical and nu- 
merical means. As we increase the rotational frequency, 
the system favors lower LLs with higher angular mo- 
menta, while the increase of SO coupling strength would 
drive the system to higher LLs with higher angular mo- 
mentum. The two trends result from the competition 
between SO coupling strength and the energy spacings 
in the LL system. We also found the winding number 
difference by 1 between the two components indicating 
the formation of half-quantum vortex. 

Weak interactions modify the single particle ground 
state dramatically. At strong SO coupling, the ring- 



structured densities develop multiple domains whose 
number shows step behavior as we increase the rotational 
frequency. At strong rotation, increasing SO coupling 
strength would favor triangular vortex lattice, deplete 
densities from the trap center, and weaken the Skyrmion 
crystal order. In the trap center the giant vortex is 
formed when SO coupling and rotation are both strong. 
Very recently in Ref. [13] we noticed the discussion of 
possible schemes for rotating SO coupled BEG system, 
exhibiting the challenge of experimental realization. Our 
findings are expected to build the understanding of the 
interplay between SO coupling and rotation. 
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